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We study the simple Hainiltonian, H = ~K{S\^ + Siz) + A5i • S2, of two, large, coupled spins 
which are taken equal, each of total spin s. The exact ground state of this simple Hamiltonian is 
not known for an antiferromagnetic coupling, A > 0. In the absence of the exchange interaction, the 
ground state is four fold degenerate, corresponding to the states where the individual spins are in 
their highest weight or lowest weight states, |t, t)5 14-! 4-) 1 Iti4-)i IJ-it)) in obvious notation. The first 
two remain exact eigenstates of the full Hamiltonian. However, we show the that the two states 
Itii)) lilt) organize themselves into the combinations |±) = -^^dt, 4,) ± |4-t))) up to perturbative 
corrections. For A > 0, the ground state is non-degenerate, and we find the interesting result that for 
integer spins the ground state is |+), and the first excited state is the anti-symmetric combination 
|— ) while for half odd integer spin, these roles are exactly reversed. The energy splitting however, is 
proportional to A^*, as expected by perturbation theory to the 2s*'' order. We obtain these results 
through the spin coherent state path integral. 

PACS numbers: 73.40.Gk,75.45.+j,75.50.Ee,75.50.Gg,75.50.Xx,75.75.Jn 



Introduction- We study the case of two large, coupled, 
quantum spins in the presence of a large, simple, easy 
axis anisotropy, interacting with each other through a 
standard spin-spin exchange coupling, corresponding to 
the Hamiltonian 



H = -K{Sl, + Si,) + \Si ■ §2. 



(1) 



We will consider K > and specialize to the case of 
equal spins Si = S2 = S. A > gives anti-ferromagnetic 
coupling while A < is ferromagnetic. The spins Si 
could correspond to quantum spins of macroscopic multi- 
atomic molecules [lH3|, or the quantum spins a macro- 
scopic ferromagnetic grains [j|, or the average spin of 
each of the two staggered Neel sub-lattices in a quantum 
anti-ferromagnet [3, l5|. 

The non- interacting system is defined by A = 0, here 
the spin eigenstates oi Siz, notationally |s, siz)(g)\s, S2z) = 
\siz,S2z), are obviously exact eigenstates. The ground 
state is four-fold degenerate, corresponding to the states 
\s, s), I — s, — s), |s, —s) and | — s, s), which we will write 
as It, t> 1 II. I) , It4) , II, t> , each with energy E = -2Ks'^. 
The first excited state, which is 8 fold degenerate, is split 
from the ground state by energy /S.E = K{2s — 1). 

In the weak coupling limit, A/A' -^ 0, it is an interest- 
ing question to ask what is the ground state and the first 
few excited states of the system for large spin S. Sur- 
prisingly, this is a yet in general unsolved problem. For 
spin 1/2, the exact eigenstates are trivially found, for 
spin 1, the problem is a 9 x 9 matrix, which again can be 
diagonalized, but soon the problem becomes intractable. 
In principle we must diagonalize a (2s -I- 1)^ x (2s -I- 1)^ 
matrix, that though is rather sparse, is not amenable 
to an exact diagonalization. For for weak coupling the 
anisotropic potential continues to align or anti-align the 



spins along the z axis in the ground state. 

As the non-interacting ground state is four fold degen- 
erate, in first order degenerate perturbation theory, we 
should diagonalize the exchange interaction in the de- 
generate subspace. However, it turns out to be already 
diagonal in that subspace. The exact eigenvalue is ±As^ 
for the two states |t,t)j 14-, 4-) ^nd then the full Hamilto- 
nian, Eqn.([T]), has total energy eigenvalue (— 2A'± A)s^. 
The full Hamiltonian can be alternatively written as 



H = ~K{Siz + Siz)±^ SizS2 



1 



(^1+^2 



■SiS+) 



(2) 
where S*,- = Six i iSiy. S^ act as raising and lowering 
operators for Siz- Thus they annihilate the states \^,^ 
),|4,,|) respectively. For the two states ItjilJiit) the 
perturbation continues to be diagonal but only in the 
degenerate subspace, with eigenvalue =fAs^ yielding, in 
first order degenerate perturbation theory, a total energy 
of (-2A'tA)s2. 

However, we will find that the states |±) = -7^(1 t 
, 1) ± I It) ) ^I's the appropriate linear combinations im- 
plied by degenerate perturbation theory. The perturbing 
Hamiltonian links the state | ± s, ^s) only to the state 
I ± s =F 1, =Fs ±1). To reach the state | =F s, ±s) from the 
state I ± s, =ps) requires one to go to 25**^ order in per- 
turbation, and s is assumed to be large. Indeed, we find 
our results via macroscopic quantum tunnelling using the 
spin coherent state path integral. Using the path inte- 
gral to determine large orders in perturbation theory has 
already been studied in field theory [6|. 

Spin coherent state path integral - The quantum (large) 
spin systems can be described by the spin coherent state 



path integral |10l . Ill | . 






(3) 



5^ is the Euchdeaii action which corresponding to dy- 
namics of particles moving on a two sphere, and which 
contains first order kinetic term, the Wess-Zumino- 
Novikov-Witten (WZNW) term for the spin degree of 
freedom [7|. 

Our two spin system, in real time, is governed by an 
action S ~ J dtC where. 



(4) 



C = J dxa.Si ■ {dtSi X d^Si) - Vi{Si) + 
dxa2S2 ■ {dtS2 X 8,82) - ^2(^2) - A^i • ^2 



where now Si — (sinfl^ cos(/)i,sin0i sin(/>i,cos0i) , i = 
1,2 are three vectors of unit norm, representing semi- 
classically the quantum spin [Jl and Ui are the values 
of each spin. In terms of spherical coordinates the La- 
grangian takes the form 



C = — (7i0i(l — cos^i) — T4 

- 0-202(1 

— A (sin 6*1 sin 82 cos(0i 



ai0i(l-cos0i)--t^i(0i,0i) 

COS02)-"^^2(e2,02) 

- (/)2) +COS0iCOS^2) • (5) 



The coefficients of the WZNW terms cti and (T2 are quan- 
tized to be a half integer, at = Ni/2. 

We consider the special case of equal spins, with ai = 
a2 = s. Our analysis is valid if we restrict our attention 
to an external potential with easy-axis, azimuthal sym- 
metry, with a reflection symmetry (along the azimuthal 
axis), as in Q, l/,(0„</.O = V{0i) = V{7r - 0,),i = 1,2. 
The potential is further assumed to have a minimum at 
the north pole and the south pole, at di = 0, and n. In 
our case the potential is explicitly 



ViS,) = V{9,,(l,,)=Ksm^d, 



(6) 



It was shown in Ref. [8|, for uncoupled spins, that quan- 
tum tunnelling between the spin up and down states of 
each spin separately is actually absent because of con- 
servation of the z component of each spin. With the 
exchange interaction only the total z component is con- 
served allowing transitions |t,i) < — >■ |i,t)- In general 
tunnelling exists if there is an equipotential path that 
links the beginning and end points. We will see that 
such an equipotential path exists, but through complex 
values of the phase space variables. 

We must find the critical points of the Euclidean action 
with t — >■ —ir, see Ref. [9|, which gives 



Ce = is(j)i{l — cos 6 

+X (sin 9i sin 6*2 cos 



The solutions must start at {6i,(j>i) — (0,0) and 
(6*2, ^2) = (tt, 0), say, and evolve to (0i,0i) = (7r,0) and 
(6»2,02) = (0,0). In Euclidean time, the WZNW term 
has become imaginary and the equations of motion in 
general only have solutions for complexified field config- 
urations. Varying with respect to (jji gives equations that 
correspond to conservation of angular momentum 



is—- (1 — cos 6*1) -I- A sin 01 sin 02 sin ((/)i 

UT 



2) = (8) 



is — (1 — cos 6*2) — A sin 01 sin 02 sin (01 — ^2) = (9) 
dr 



While varying with respect to 0; gives the equations: 



901 ^ ^ 002 
Adding Eqn's ^ and ^ we simply get 
d 



dr 



(cos 01 + cos 02) = 0. 



Hence 



cos 01 -f- cos 62 = I = 



(10) 



(11) 



(12) 



where the constant I is chosen to be zero using the initial 
condition 0i = 0, 02 = tt. We can now eliminate 02 from 
the equations of motion and writing = 0i , = 01 — 02 
and $ = 01 -H 02 and taking Vi{e,) = V{e,) = l/(7r - 0,) 
we get the effective Lagrangian: 



£ = is$ - is0 cos + C/(0, 0) 



where 



U (0, 0) = 2V (0) -I- A (sin^ cos - 



(13) 



+ A (14) 



is the effective potential energy, where we have removed 
a constant A so that the potential is normalized to zero at 
= 0. The first term in the Lagrangian is a total deriva- 
tive and drops out. The equations of motion become: 



isf smf 



9^(0,0) 
0(1) 



(15) 
(16) 



These equations have no solutions on the space of real 
functions 0(t), 0(t) due to the explicit i on the left hand 
side. The analog of conservation of energy follows im- 
mediately from these equations, this is easily derived by 
multiplying (fT5|) by 9 and ([T5)) by and subtracting, this 



gives: 



V{9i) + isMl - cos 02) + V{92) 

'1 - 02) + cos 01 cos 02) . (7) 



dU (0, 0) 
d^ 



i.e, t/ (0, 0) = const. = (17) 



The constant has been set to using the initial condition 
9 = 0. From Eqn. p4| . we have 



(18) 



U {d, (j)) = (2A' + A (cos + 1)) sin^ 9 = 



implying {2K + A (cos + 1)) = as sin^ 9 ^0^ which is 
required for a non-trivial solution. Then 



C0S( 



2K 



1 



(19) 



and we see that (j) must be a constant. This is not valid in 
general, it is due to the specific choice of the external po- 
tential Eqn. ^. Since K > |A| we get | cos0| > 1, which 
of course has no solution for real (p. We take cf) = (l)ji+i(pj 
which gives cos (f> = cos (f>R cosh (pi — i sin (p^i sinh (pi . As 
the RHS of Eqn. (fT9|) is real, we must have either 0/ = 
or (f)ft = nil or both. Clearly the (pi = cannot yield 
a solution for Eqn. (|19p . hence we must have 0^ = mi. 
As we must impose 27r periodicity on pu only ri = or 1 
exist. Then we get 

r-(^ + i) ifA>o 

cos0H-ircosh0.^|/^^_/^ ifA<0 ^'^ 

Thus ?7 = 1 for A > and n = for A < allowing for 
the unified expression 

cosh0, = ?^. (21) 

Eqn. (fT6|) simplifies to 



is- 



sm( 



— Asin(/) = — iA(— l)"sinh0/ = z|A| sinh0/ 

(22) 

as A(— 1)" = — |A|. Eqn. ((2T|) has two solutions positive 
(pi corresponds to the instanton, (^ > 0), and negative pi 
corresponds to the anti-instanton, {9 < 0). The equation 
is trivially integrated with solution 



/ ^(r-ro)A 



2arctan(e"(^"^°^ 



(23) 



where lj = (|A|/.s) sinh0/ and at r = tq we have 9{t) = 
tt/2. Thus 9{t) interpolates from to tt as r = — oo — > oo 
for the instanton and from tt to for an anti-instanton. 
Using (p = and (fT7|) we see that the action for this 
instanton trajectory, let us call it 5*0, simply vanishes 



Sn = 



drC = 0. 



(24) 



So where could the amplitude come from? We have 
not taken into account the fact that (j) must be trans- 
lated from (j) = (any initial point will do, as long as 
it is consistently used to compute the full amplitude) to 
(j) = nn + icpi before the instanton can occur and then 
back to (j> = after the instanton has occurred. Nor- 
mally such a translation has no effect, either the change 



at the beginning cancels that at the end, or if the action 
is second order in time derivative, moving adiabatically 
gives no contribution. But in the present case, before the 
instanton occurs, 9 = 0, but after it has occurred, 9 = tt. 
As p is multiplied by cos 9 in the action, the two contri- 
butions actually add, there is a net contribution to the 
action. Indeed the change of the full action for the com- 
bination of the instanton and the changes in (j) is given 
by 



f>n-K-\-i(f>i f-O 

AS ~ / —isdcj) COS 6\0^Q -\- So -\- / —isd(j)cos( 

Jo J n7T-\-i(f)j 



-is2n7r + 2s(j)j 



(25) 



We will use this information to compute the following 
matrix element, using the spin coherent states \9, 0) and 
the lowest two energy eigenstates \Eo) and \Ei): 

{9f,<j>f\e-P"\9,,cl>,) = e-^''°{9f,<j>f\Eo){Eo\9,,<j>,) 

+e-^'''{9f,<l>f\Ei){Ei\9,,<P,) + --- (26) 

On the other hand, the matrix element is given by the 
spin coherent state path integral 



{9f,^f\e-^"\9,,^,)=Af V9VPe 



0f-4'f 



(27) 



The integration is done in the saddle point approxima- 
tion. With {9i,pi) = (0,0) corresponding to the state 
|t,i) aud {9f,<f)f) = (tt, 0) corresponding to the state 
I ^, t)we get, with a mild abuse of notation 



(;,t 



.-PHy 



|t,;)=AAe-^^K/?(l + ---) (28) 



where k is the ratio of the square root of the determi- 
nant of the operator governing the second order fluctu- 
ations, except for the time translation zero mode, about 
the instanton and that of the free determinant. It can in 
principle be calculated, but we have not done this. The 
zero mode is taken into account by integrating over the 
position of the occurrence of the instanton giving rise to 
the factor of /3. A/" is the overall normalisation including 
the square root of the free determinant which is given 
by Ne~^°^ where Eq is the unperturbed ground state 
energy and A'^ is a constant from the ground state wave 
function. The result exponentiates, but since we must 
sum over all sequences of one instanton followed by any 
number of anti-instanton instanton pairs the total num- 
ber of instantons and anti-instantons is odd, and we get 



'k;3 ^ sinh {c-^^kP) 



(29) 



Given AS* = —is2mT + 2s(pi and solving Eqn. (PTjl for 

4)1 for a:> iiai 



arccosh 



W)-( 



AK 

'^'w 



(30) 



gives |1^ 



-,-AS 



^is2TZ—2s(f>i 



if A>0 






2s 

if se Z 



- m'^'^ 



4K 



2s 



if A<0 



Then we get 



a,t 



= -/3ff 



it,;) = ± 



i /-Mr 

2 



'k/3 



1 _(-Mf 
— e V4K/ 



'kP 



-1/2 
(31) 

(32) 
where the — sign only apphes for the case of anti- 
ferromagnetic couphng with half odd integer spin, A > 
0, s = Z + 1/2. An essentially identical analysis yields 

a,tie-^^i;,t> = (t,;ie-^^it,;> 

= fle{m"'-f^ + le-i^f-ANe-^''°. (33) 

These calculated matrix elements should now be com- 
pared with what is expected for the exact theory: 

(34) 
and say 

(;,t|e-'^^|;,t) = e-f'^^-~i^^^{U\Eo){Eo\iA) 

+ e-P^''°+-^^^^{U\Ei){E,\U) 

(35) 

The energy splitting can be read off from this result 
AE = Ei-E2 = 2 ( i^ j K (36) 

and our main result follows, the low energy eigenstates 
are given by 

|£^o) = ^(ii>t> + it,;» \Ei) = ^(i;,t> - It,;)) 

(37) 
for A < (although here the energy eigenstates are 
should be lE^) and I-E4)) and A > for s G Z, while 
for the anti-ferromagnetic A > case with half odd inte- 
ger spin we get 



\Eo) = -^(i;,t) - it,;» \E^) = -^(i;,t> + it,;». 

(38) 



This understanding of the ground state in the anti- 
ferromagnetic case is our main result. 

Conclusions- We have found the low energy eigenval- 
ues and the corresponding eigenstates for the Hamilto- 
nian of two equal, large, spins interacting with an easy 
axis anisotropy and a standard exchange interaction, the 
latter which is considered as a perturbation. We find that 
the two states |4,,t), It,!) reorganize into the symmetric 
and the anti-symmetric superposition because of quan- 
tum tunnelling transitions. These transitions correspond 
to the 2s'^ order effects in perturbation theory. The sym- 
metric combination is the lower energy state for integer 
spin while the anti-symmetric state is the the lower en- 
ergy state for half odd integer spins. These states are 
respectively the ground states for an anti-ferromagnetic 
coupling. 
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